FIRST DERIVATIVES ESTIMATES FOR 
FINITE-DIFFERENCE SCHEMES 



ISTVAN GYONGY AND NICOLAI KRYLOV 

Abstract. We give sufficient conditions under which solutions of 
discretized in space second-order parabolic and elliptic equations, 
perhaps degenerate, admit estimates of the first derivatives in the 
space variables independent of the mesh size. 



1. Introduction 

This is the first article out of a series of two devoted to estimating 
space derivatives of solutions of discretized in space second-order par- 
abolic and elliptic equations. We allow equations to degenerate and to 
become just first-order equations. In the present article we only deal 
with the first-order derivatives. In the second part of this project we 
will prove higher-order derivatives estimates and apply them to show- 
ing a method of accelerating finite-difference approximations to any 
given rate for equations in the whole space. 

Numerical approximations for linear and quasilinear partial differ- 
ential equations is a rather old and well developed area. We refer the 
reader to [2] and the references therein, following which one can track 
down original papers by D. Aronson, L. Bers, R. Courant-K. Friedrichs- 
H. Lewy, J. Douglas, F. John, O. Ladyzhenskaya, P. Lax, H. Levy, L. 
Liusternik, I. Petrovskii, and many many others to which we only add 
[8] and two more papers |I] and [9] where discrete methods are applied 
to approximate stochastic partial differential equations. 

A major difference of this article from all above mentioned ones is 
that we focus on investigating the smoothness of approximating solu- 
tions rather than on convergence only. For each point x G M d we move 
the original grid in such a way that x becomes a grid point. This allows 
us to define the approximate solution in all of M. d rather than only on 
the grid and we investigate how smooth the approximating solution 
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is with respect to x. We estimate true derivatives rather than their 
difference approximations. 

Estimating the sup norms of the first-order derivatives for solutions 
of finite-difference schemes for linear and fully nonlinear second-order 
degenerate equations plays a major role in estimating the rate of con- 
vergence of approximating solutions to the true solution in the sup 
norm. The most general results for fully nonlinear equations concern- 
ing the rates can be found in [5] and in the references therein. A re- 
cent development in the issue of estimating the Lipschitz constant and 
second-order differences for approximating solutions for fully nonlinear 
equations without applications to estimating the rate of convergence 
is presented in [6]. Before, the Lipschitz constants and higher order 
derivatives estimates were obtained in [2] for time-space discretization 
of linear degenerate parabolic equations and in [1] for fully nonlinear 
equations. They are also applied to estimating the rate of conver- 
gence. In a sense the present article is close to [2]. However, here we 
only deal with the first-order derivatives estimates and for equations 
discretized only in the space variable. We introduce a new type of suf- 
ficient conditions for obtaining the estimates (see Assumption 12.31 and 
the discussion in Section These conditions are much weaker and 
more detailed than the corresponding ones in [2]. Our method is also 
somewhat different. Instead of considering just the sum of squares of 
the difference increments along the mesh we add to it the square of the 
full gradient with a small constant factor. This allows us to estimate 
the gradient. 

In this connection it is worth noting that such an estimate for finite- 
differences approximations of the first-order directional derivatives in 
x is claimed in Theorem 4.1 of [2] under some conditions, which are 
always satisfied if the equation is uniformly nondegenerate even if c 
(see (12.21) ) is small. However, in this case, actually, the result of The- 
orem 4.1 of [2] is only proved for the derivatives along the mesh. This 
is rather harmless if the vectors on the mesh span the whole space, 
but excludes the cases when the mesh lies in a subspace, which hap- 
pens, for instance, if we are dealing with, say uniformly nondegenerate 
equations whose coefficients depend on a parameter and we want to 
estimate the the finite-differences of their solutions with respect to the 
parameter by considering it as just another space variable. In that case 
no second-order derivatives with respect to the parameter enters the 
limit equation, the assumption that it is uniformly nondegenerate with 
respect to the original space variables does not help, and we need to 
have c be large in order to rely on Theorem 4.1 of [2]. Our results are 
free from this flaw, see Remarks 15.31 and 16.51 
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To understand faster the method and the results of the article, we ad- 
vise the reader concentrate only on the parabolic case and assume that 
the limit equation is uniformly nondegenerate. Then apart from the 
Lipschitz continuity nothing else (see Remarks 12.11 and 12. 3j) is required 
for Theorem 12.11 to hold. 

For the general equations our conditions (see Remarks 16.21 and 16.61) 
capture the main features of the corresponding conditions known from 
the theory of PDE. Namely, roughly speaking, we need the first-order 
derivatives in any direction along the mesh of the coefficients to be 
dominated either by the diffusion coefficients along the same direction 
or by c or else by the drift term if it is sufficiently "monotone" . It is 
worth noting that along the way we discover the necessity of the dif- 
fusion coefficients to have a special form and the usefulness of adding 
a diffusion term with a coefficient proportional to the mesh step into 
approximating equation. This reminds the method of artificial diffu- 
sion, although, as far as we understand, the artificial diffusion is usually 
added to the original differential equation. 

Of course, in the same way as in [2], the results of the present article 
lead to the rate of convergence of order h l l 2 of approximating solutions 
to the true solution. However, for brevity we do not say more about 
this issue only adding that in general our finite-difference equations 
need not be related in any way to a partial differential equation. 

Our main results are collected in Section [21 which also contains the 
proofs of all of them but Theorem 12.11 which is proved in Section 
[31 In Section [2] we also state, in a special case, without proof one of 
the main results of the continuation of the present paper. Section [4] 
contains a discussion of our assumptions concerning the structure of the 
finite-difference equations under consideration. The point is that our 
equations do not contain mixed second-order differences and in Section 
H]we explain that this is "almost" the most interesting case. The final 
Sections [5] and are devoted to a rather long and detailed discussion 
of the somewhat formally stated Assumption 12.31 and showing that it 
is natural in many cases alluded to above. 



We take some numbers h,T £ (0, oo) and in a cylindrical domain 
consider the integral equation 



Jo 

for u, where g(x) and f(s,x) are given real-valued Borel functions of 
x = (xi,...,Xd) £ M d and (s,x) £ H T : = [0, T] x M. d , respectively, and 
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L is a linear operator given by 

Lip(t, x) = Lh(f(t, x) = L°Lp(t, x) — c(t, x)(p(x), 

L°<f(t,x) = L° h <f(t,x) = - ^2 qx(t,x)5 X (p(x) + ^2p x (t,x)5 X ip(x), 

AeAi AeAi 

(2.2) 

for functions ip on M. d . Here Ai is a finite subset of M. d such that Ai, 
px(t,x), q x (t,x) are real-valued functions of (t,x) G Ht given for each 
A G Ai, and 

5 x ip(x) = 6 ht \<p(x) = \{(f{x + hX) - <p{x)), A G Ai. 

h 

Let m > be an integer and let Ki G [1, oo) be a constant. Introduce 

Xx = Xh,x :=qx + hp x . 
We make the following assumptions. 

Assumption 2.1. The functions p, q, c, /, and g and their derivatives 
in x up to order m are bounded on Ht and continuous in x. 

Assumption 2.2. For all (t,x) G H T and A G Ai, 

XaM >0. 
There exists a constant Co > such that c> cq. 

Remark 2.1. The above assumption: c > Co > 0, is almost irrelevant 
if we only consider (12.11) on a finite time interval. Indeed, if c is just 
bounded, say |c| < C — const, by introducing a new function v(t,x) = 
u(t,x)e~ 2Ct we will have an equation for v similar to ( 12. ip with L°v — 
(c + 2C)v and fe~ 2Ct in place of Lu and /, respectively. Now for the 
new c we have c + 2C > C. 

Remark 2.2. Introduce the following symmetry condition: 
(S) We have Ai = — Ai and = Q-x on Ai. 
Obviously under condition (S) we have 

h" 1 1x(t,x)5 x tp(x) = (1/2) J2 qx(t,x)A X ip(x), 
AeAi AeAi 

where 

A A (/?(x) = h~ 2 (Lp(x + hX) - 2ip(x) + (p(x - hX)). 

Take a function r x defined on A x taking values in [0, 1], and for A G Ai 
introduce the operators 

T\<p(x) = T KX ip{x) = tp(x + hX), 5 X = r x h~ l (T x - 1). 
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It is worth noticing that in most applications we take T\ = 1 on Ai. 
However, there are cases (see Remark 15.31) in which it is useful to have 
some flexibility in changing T\. 

For uniformity of notation we also introduce A 2 as the set of fixed 
distinct vectors i 1 , ...,£ d none of which is in Ai and define 

8# = hj* = T Di := r d/dxi, T(i = 1, A = Ai U A 2 , 

where To G [0, 1] is a fixed constant. Observe that we allow r to be zero 
in order to cover some results from [2]. 
For ji G A we set 

Qif = h7 l ^2 <?a<W, = h~ l ^2 (<W)£\<£, 

AeAi AeAi 

P( P= Pa<W, Pp<p = ^ (^Pa)*aV, 
AeAi AeAi 

Below B(M. d ) is the set of bounded Borel functions on M. d and ^ 
is the set of bounded operators K, = JC(t) mapping B(R d ) into itself 
preserving the cone of nonnegative functions and satisfying /CI < 1. 
We will often make use of the simple fact that for any /Q, /C 2 G ^ and 
nonnegative functions a, (3 on R d , 

a/Ci + j3JC 2 = (a + /3)/C 3 

with 

^ := ^/Ci + ^/C 2 G S, (g:=0). 

Assumption 2.3. We have m > 1 and there exist a constant 5 G (0, 1] 
and an operator K, — ]Ch G .ft, such that 

2 ^(A»L° A T A <p < ^ Q(A» + K x Q(ip) + 2(1 - 5)c/C( ^ l^| 2 ) 
AeA AeA AeA 

(2.3) 

on Ht for all smooth functions if, where 

Q(<p) = E ^IM 2 - 

z^eAi 

It is worth noting that Assumption 12.31 is automatically satisfied if 
qx and p\ are independent of x. There are a few more cases when it is 
satisfied as well. We discuss some of them here and in Section [5] only 
mentioning right away three situations. 
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Remark 2.3. Let Assumptions I2.ll and 12.21 hold with m > 1. Assume 
that Ai = — Ai, Dq_\ — Dq\ and q\ > K for all A G A 1; where n > 
is some constant. Then Assumption 12.31 is satisfied as well with 5 as 
close to 1 as we wish, with t\ = 1 on Ai, appropriate r > 0, K\, unit 
JC, and all small h. 

We will prove this remark at the end of this section. In Remark 16.41 
we show that if we have m > 2 and the symmetry condition (S) is 
satisfied, then in the above remark the condition k > can be replaced 
with k = 0, provided that cq is large enough (this time we need not 
assume that h is small). In Remark 16.41 we also show that the condition 
m > 2 can be replaced with m = 1 provided that ^fq~\ are Lipschitz 
continuous in x with a constant independent of t. In that case again 
Assumption 12.31 is satisfied for 5 = 1/10 and appropriate K\,tq > 0, 
provided that c is large enough. 

As we have seen in Remark 12.11 the condition that Co be large is, 
actually, harmless as long as we are concerned with equations on a 
finite time interval. 

Fix a domain Q C M d and introduce 

Q° = { x e Q : x + Xh G Q VA G AJ, 5Q = Q\Q°, 
Q T = [0, T] x Q, Q° T = [0, T] x Q°, 6 X Q T = [0, T] x 6Q, 
S'Qt = ({0} xQ)U 5 X Q T . 
Our first main result is formulated as follows, where by Du we mean 
the gradient of u with respect to x. Observe that the main case that 
Q = M. d is not excluded and in this case assumption (ii) of Theorem 
12.11 below can be checked on the basis of Theorem 12.31 A typical and 
the most reasonable application of Theorem 12.11 when Q is a proper 
domain is the case that r = 0. 

Theorem 2.1. (i) Let Assumptions \2.1\ throuah \2.3\ be satisfied and let 
u be a bounded function on Ht satisfying (12.11) in Qt- 

(ii) Assume that u and Du are bounded and continuous in Qt- 

Then in Qt we have 

\u\ + t \Du\ + U < N(F 1 + sup(|u| + t \Du\ + £/)), (2.4) 

S'Qt 

where 

U=(^\Sxu\ 2 ) 1/2 , F 1 = S a V {\f\ + \Df\), 

and N depends only on S,cq,Ki, sup^ T \Dc\, and 

|Ai| 2 := E l^l 2 . 
AeAi 
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Remark 2.4. We will see from the proof that, if To = 0, then Theorem 
12.11 holds without the assumption that Du exists let alone continuous. 

In case Q = M. d assumption (ii) of the following result is often satis- 
fied due to Theorem 12.31 

It is worth noting that if Q = ~R d , then 5Q = and for any function 

<p we set 

sup ip := 0. 



Theorem 2.2. (i) Let Assumptions \2. 1\ through \2. 3\ be satisfied. Sup- 
pose that q\, p\, c, and f are independent oft. 

(ii) Assume that in R d there exists a bounded function u = u(x) 
which is bounded and continuous in Q along with Du and such that 

Lu + / = in Q. 

Then in Q we have 

\u\ + t \Du\ + U < N(F 1 + supO| + t \Du\ + U), 

SQ 

where U is the same as in Theorem \2.1\ 

F 1 = BU P (|/| + |D/|) J 

and N depends only on 5,c ,Ki, sup M d \Dc\, and |Ai|. 

Proof. Take v = cq/2, so that c — v > cq/2, and observe that in Qt 
the function v(t,x) := u{x)e vt satisfies 

^-v = L°v - (c - u)v + e ut f. (2.5) 

By Theorem 12.11 for x G Q and obvious meaning of V we have 
e vT {\u{x)\+r \Du{x)\ + U{x)) = \v(T,x)\+r \Dv(T,x)\ + V(T,x) 

< Are^[Fi + sup(|u|+T |£>u| + [/)] + iVsup(M + t \Dv\ + V)(0, y). 

SQ Q 

By multiplying the extreme terms by e~ uT and letting T — > oo, we get 
the result. The theorem is proved. 

Remark 2.5. It is worth noticing that in the above theorems it suffices 
that (12.31) be satisfied only in Qj,. 

Theorem 2.3. (i) Let Assumption \2. II be satisfied. Then there exists 
a unique bounded solution u of (12. ip in Ht- Moreover, all derivatives 
in x of u of order < m are bounded and continuous in Ht- 
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(ii) Let Assumptions \2. il through \2.3\ be satisfied. Suppose that q\, 
Px, c, and f are independent oft. Then there exists a unique bounded 
solution u = u(x) of the equation 

Lu + f = in R d . (2.6) 

Moreover, u and Du are bounded and continuous in W 1 . 

Proof, (i) Let C m be the space of functions on M. d which are bounded 
and continuous along with all derivatives up to order m. We endow 
C m with an appropriate sup norm and in so obtained Banach space, 
denoted again by C m , consider the equation 

u(t) = g+ [ (A{s)u(s) + f{s))ds } 



where f(s) — f(s, ■) and A(s) are operators in C m given by 

A(s)ip(x) = h" 1 ^ X\(s,x)5 x ^(x) - c(s,x)ip(x). 
AeAi 

Owing to Assumption 12.11 

||A(s)</?||c- < Ny\\ C m 

with N independent of s and </?. Hence, our result is a direct conse- 
quence of the general theorem about ordinary differential equations in 
Banach spaces (for proving uniqueness we take m = 0). 

(ii) By assertion (i) for any T there exists a unique bounded and 
continuous in Ht solution v(t, x) of the problem 

d 

—v(t,x) = (L + u)v(t,x) t>0, v(0,x) = f(x), 

where v = Cq/2. In addition, Dv is bounded and continuous in Ht for 
each T. By Theorem 12.11 v and Dv are bounded and continuous in 
Hoq. Define 

POO 

u(x) = / e~ ut v(t,x)dt. 
Jo 

Then the rules of differentiating under the integral sign and the dom- 
inated convergence theorem show that u and Du are bounded and 
continuous. Furthermore, integrating by parts, we see that 

Lu(x)= I e~ ut Lv(t,x)dt = / e~ ut [—v{t,x)-uv{t,x)]dt = -f(x), 
Jo Jo fit 

so that u satisfies ( 12. 6p . 

To prove uniqueness of bounded solutions of (12.61) we use Lemma 
13.11 which is proved in Section [31 If w is the difference of two bounded 
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solutions of (12.61) . then v(t, x) := w(x)e cot satisfies ( 12. 5p with v = Co 
and / = 0. Since c — v > 0, by Lemma [3. II 

sc) < supw + (0, = supw + , w(x) < e~ c °*supu; + 

and by letting t — > oo we obtain w < 0. The same inequality holds for 
— w, so that to = 0, which proves uniqueness and finishes the proof of 
the theorem. 

Remark 2.6. Let Assumption 12.11 hold. Then it is easy to see that for 
the bounded solution u of (12.11) in Ht for each fixed hi > 

sup \u(t,x)\<N sup \f(t,x)\, 
(t,x)eH T (t,x)eH T 

holds for all h > hi, where is a constant independent of h. 

Remark 2.7. A simple inspection of their proof shows that the above 
theorems remain valid if /C( XIaga l^^l 2 ) m Assumption 12. 31 is replaced 
by 

sup |<W| 2 - 

X AGA 

We do not know how much can be gained by such weakening of As- 
sumption 12.31 On the other hand, in a subsequent article we will see 
an advantage of using operators JC e 

Now we state without proof one of the main results of the forth- 
coming paper [5j. As we know from Remark 12.31 and Theorem 12.31 (ii). 
under the conditions of Theorem 12.41 (see below), for each h > 0, there 
exists a unique bounded solution Uh of 

L h u + / = in R d . 

For a fixed integer k > set 

k 

where 

(b ,bi,...,b k ) := (1,0,0,... ,0)^ 
and V~ x is the inverse of the Vandermonde matrix with entries 
yij .- 2 -(i-l)(i-l) ) ij = 1, ...,&+!. 

Consider also the equation 

Cv + f = in R d (2.7) 

with 

£ := OijDiDj + biDi - c, (2.8) 
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aij(x) := (1/2) ^ Qx^XiXj, h := ^ Px(x)K- 

AgAi AgAi 

One of our main theorems from [5] in a special case reads as follows. 

Theorem 2.4. Let m > 3(k + 1) for some integer k > 0. Let Assump- 
tions {KJl and the symmetry assumption (S) be satisfied. Also 
assume that q\(x) > k for all A G Ai, i G l rf for some constant n > 0. 
Then there is a unique bounded solution v to (12.71) and 

\v h (x)-v(x)\ < Nh k 

for all x G M d , h G (0, ho], and every h > 0, where N is a constant 
depending only on ho, k, m, k, Co, |Ai|, and on the sup norms of the 
derivatives of q\, p\, c, and f up to order m. 

We obtain this result in [5] by showing that the derivatives of v h in 
h up to order k + 1 are bounded functions of h G (0, ho], which we will 
prove via our estimates on the derivatives of Vh in x. The reader may 
wonder why we do not choose the straightforward way of estimating 
the derivatives of Vh in h via an 'explicit' formula for v. To test this 
approach we suggest the reader try to estimate dvh/dh directly for 

oo n 

v h (x) = ^f(x) + £ h\^T +1 Ef(x + £ he*), 

n=l i=l 

where are independent random variables taking 1 and —1 with prob- 
ability 1/2, without noticing that Vh is the bounded solution of 

-jj>(u(x + h) - 2u(x) + u(x - h)) - u(x) + f(x) = 0, x G R. 

We finish the section by proving the assertion in Remark 12.31 Clearly, 

2j2(hv)L° x T xV = I 1 + I 2 , 

aga 

with 

Ji := 2 (h<p)L° x T x <p, h:=2j2 (M^aV- 

AgAi AgA 2 

We take t x = 1 on Ai and notice that due to the symmetry of Ai and 
the symmetry of Dq x in A 

h = 2 ( S w) L °x<P + 2hJ2 (Sx<p)L° x 5 x cp 

AgAi AgAi 

= (<W)(5a<?m) a ^ + 2 Y (h^iSxP^S^ 

A,/xgAi A,^gAi 

+2 M(SxX^S X ip =: + lf ] + J{ 3) , 

A,/iGAi 
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d 



3=1 MSAi 




2r 2 ^^(D jf )(%)^ 



A*eAi 



By Young's inequality, taking into account that xa > K /2 > for 
sufficiently small h, and that c > Co > 0, we have 



A j) < (1/3) ^ Q(5 x <p) + JVQ(^) for j = 1, 3, l[ 2) < NQ(ip), 



where N is a constant depending only on k, the number of elements 
in Ai and on the supremum norm of the gradients of p\ and q\ in x. 
Summing up these inequalities and taking r > sufficiently small we 
get (12.31) with K\ = 2N, unit operator /C, and with 5 as close to 1 as 
we wish. 



If Q = M. d , Theorem 12.31 (i) shows that equation (12. ip has a unique 
bounded continuous solution u for which the partial derivatives in x G 
M. d up to order m are bounded continuous functions of (t, x). However, 
the bounds, which can be extracted from the proof of Theorem 12.31 
for these derivatives depend on the parameters h and T. Our aim is 
to show the existence of bounds, independent of h and T if m = 1 
and in addition to Assumption 12.11 Assumptions 12.21 and 12.31 also hold. 
We will obtain such estimates by making use of the following version 
of the maximum principle. It is probably worth drawing the reader's 
attention to the fact that the assumption that c has certain sign is not 
used in Lemma 13.11 

Lemma 3.1. Let Assumption \2. 1\ with m = be satisfied and let x\ > 
for all A 6 Ai. Let v be a bounded function on Qt, such that v(-,x) is 
measurable for any x G Q and the partial derivative D t v := dv(t,x)/dt 
exists in Q^. Let F(t) > be an integrable function on [0,T], and let 



AGAi 



IP < (1/3) £ Q(M + rlNc^c I'M 2 , 




4 2) <T 2 Nc^cJ2\5xV\ 2 , 



3. Proof of Theorem 12.11 
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C(t,x) > be a bounded function. Assume that for all (t,x) G Q° T we 
have 

D t v < Lv + Cv + + F, (3.1) 

where v(t) = swp{v(t,x) : x G Q}. 
Then in [0, T] we have 

v(t) < G(t)e ut + f F{s)e u(t - s) ds, (3.2) 
Jo 

where 

^ := sup(C — c), G(t) = sup e~" s v + (s, y), 
Q° T (s,y)e8'Q t 

s'Qt = ({o} x q) u 5 x g t . 

Proof. First assume that C = c = F = G = 0. In that case introduce 
v(t,x) = v(t,x)e Nt , where 

N = N h = sup h~ l ^ xx- 

Qt AeAi 

Observe that w(t,x) = v(t,x)e~ Nt , 

e- m {D t v{t,x) -Nv(t,x)) < e- m Lv(t,x), 
D t v(t,x) < (L + N)v(t,x) 

in Q T , and in 



v(t,x) < Jv(t,x) := J Q o / (L + iV)(/ QT i;)(s,a;)ds. 

./o 

Obviously, due to the choice of iV and the assumption that Xa > 0, if 
u\ > «2 o n Qt ; then Jw^ > Ju 2 on Qt- It follows that, < J k v for 
any fc. 

Observe that for any bounded function / 

sup\Jf(t,-)\<N' [ sup\f(s,-)\ds, 
Q Jo Q 

where N' is independent of / and t. It follows easily that for any 
bounded function / we have J k f — > uniformly on Qt- Hence v < 
and v < in Qt- 

Now we consider another particular case in which G = c = and 
C = const > 0. Then observe that the function 

w{t,x) = [ (Cv+ + F)(s)ds 



satisfies 

D t w = Lw + Cv + + F. 



FIRST DERIVATIVES ESTIMATES 13 

Therefore v := v — w satisfies D t v < Lv. In addition v < on 5'Qt- 
By the above, on Qt we have v < 0, that is 

v(t,x)< [ (Cv+ + F){s) ds v+(t)< [ (Cv+ + F)(s)ds, 
Jo Jo 

and (13.21) follows by Gronwall's inequality and the fact that v = C. 

Now we allow and variable C but still assume that G = 0. In 
that case take a large constant M so that M > c and M + v > and 
for v(t, x) = v(t, x)e Mt write 

e- Mt (D t v(t,x) - Mv(t,x)) < e- Mt Lv(t,x) + C(t,x)v+(t) + F(t). 

Dropping obvious values of arguments and introducing 

v(t) = snpv(t, x) (= ve Mt ), 
xeQ 

we find 

D t v < L% + (M — c)v + Cd + + e Mt F < L°v + (M - c)d+ 

+C% + + e Mt F < L% + (M + + e Mt F. 
It follows by the above that 

v(t)e Mt < f e Ms F(s)e {M+ ^ t ~ s Us, 
Jo 

which is equivalent to (13. 2p . By the way, notice that so far we have not 
used the fact that C > 0. 

Now comes the general case in which we set 

w(t, x) = v(t, x)e~ ut - P, (3= sup v + (s,y)e~ us . 

(s,y)£5'Q T 

Simple manipulations show that (13.11) becomes 

D t w < L°w - c(w + 13)- u(w + /3) + Ce~ ut v + + Fe~ ut , 

where 

e~ ut v + < w + + j3, w+(t) := sup w + (t,x) 

Q 

and, owing to the assumption that C > 0, the definition of u, and the 
fact that f3 > 0, 

Ce~ ut v + < Cw + + Cf3< Cw + + {c + v)f3. 

It follows that 

D t w < L°w - (c + v)w + Cw+ + Fe~ v \ 
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Since w < on 5'Qt and sup(C — (c + u)) = 0, by the above cases we 
have for t < T that 

w{t) < [ F{s)e~ us ds, 
Jo 



u(t)<e"* sup v+{s,x)e- vs + e vt F(s)e- ys ds. (3.3) 

{s,y)e8'Q T Jo 

We can put here t = T and then, by using certain freedom in choosing 
the end of the time interval, we can, actually, set T = t in (13.31) . Then 
we arrive at (13. 2p for all t < T. The lemma is proved. 

Corollary 3.2. Under the conditions of Lemma \3. II if v < 0, then 

v(t) < sup v + (s, y) + sup F. 
(s,y)eS'Q t [o,t] 

Proof of Theorem 12.11 Introduce 

V = u 2 , V 1 = ^|5 AM | 2 , V k {t) =supV k (t,x), 
AeA xeQ 

recall that Fi is introduced in the statement of Theorem 12.11 and set 

G = sup(|it| + r \Du\ + U). 

S'Qt 

By Corollary 13. 21 applied to u or — u from the assumption that c > Cq 
we obtain that 

K, = |«| 2 <(G + Co 1 F 1 ) 2 . 
Now we use the formula 

<5 A 0M = {hil>)T\<p + if>S X (p = (5 x ip)^P + #a¥> + h(5 x if})5 X (p, A e Ai, 
to get 

L°((p 2 ) = 2<pL°<p+Q(<p). 

In particular, 

LVx = 2 ^(5 A n)L°5 A n + ^ Q(5 A n). 

AeA AeA 

We observe that 

L°5 x u = 5 x L°u - L° x T x u 

if A G A, and in Q° T 

6 x L°u = D t 5 x u + c5 x u + (5 x c)T x u - 5 x f. 
Then in Q° T we find 

L°Vi - DtVr - 2cVi = 2(M 

AeA 
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+2 ^(M {{hc)T x u - 8 x f) - 2 ^(ML°T A u . 

AgA AGA 

We use Assumption 12.31 to conclude 

L°Vi - D t Vi - 2cV! > -KiQ(u) - 2(1 - S)cVi + /, (3.4) 

where 

I:=2^2(hu)({S x c)T x u-6 x f). 

aga 

Notice that by Young's inequality 

aga aga 
where, for each A G A l5 

f h 

5 x c(t, x) = /T 1 / XiDic(t, x + A0) d0, 
./o 

\5 x c(t,x)\ 2 < \r x \\ 2 h~ l / \Dc(t,x + \6)\ 2 d9 <\t x \\ 2 sup \Dc\ 2 , 
so that 

Y,(^?<^v\Dc\ 2 {l + \\ l \ 2 ) 1 

AGA ^ T 

2 ^ \{6 x u){(6 x c)T x u\ < 5 2 V 1 + NV < 5 2 V t + N(F 2 + G 2 ). 
aga 
Similarly, 

2Y,\{5xu)5 x f\<5 2 V l + NF 2 . 

aga 

Hence (J33D yields 
LVi-AK-2cV[ > -2(5 2 + c-5c)V l -K 1 Q(u)-N(F 2 + G 2 ). (3.5) 
Next, 

L°(u 2 ) - 2cu 2 - D t (u 2 ) = 2u(L°u - cu - D t u) + Q(u) 

= -2uf + Q(m) > —N(F 2 + G 2 ) + Q(u). 
This, (13.51) , and the fact that S 2 + c — <5c > show that for 

H/ := V x + ^u 2 , H/ = sapW(-,x), 

Q 

we have 

L^H/ - 2ciy - D t W > -2(5 2 + c - 8c)W - A(F 2 + G 2 ). 
Now we want to use Corollary 13.21 Set 

v := sup[2(5 2 + c - 5c) - 2c] = 25sup(5 - c) < 25(5 - c ). 
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If 

v < — Sc (3.6) 

then by Corollary 13. 21 we get W < N(F 2 + G 2 ), which obviously implies 
(12.41) . Finally, observe that, if Assumption 12.31 is satisfied with a 5 = 
5q > 0, then it is also satisfied with any 5 G (0, do] and the same K±, JC. 
Hence by modifying 5 if necessary, so that 5 < c /2, we satisfy (13.61) 
thus proving the theorem. 

The following remark will be used in a subsequent paper when we 
will be estimating higher order derivatives of u. 

Remark 3.1. Suppose that, instead of Assumption |2~UI Assumption 15. II 
(see Section [5]) is satisfied. Then a simple inspection of the above proof 
shows that in place of (I3.5P we would have 

L°Vi - D t Vi - 2cVi > -2(5 2 + c - 5c) V l 

AeA 

4. Some issues related to the convergence L h — > £ 

There is a natural question about the relation of the finite-difference 
operators L = Lh with partial-differential operators. We certainly want 
to apply the results of the present article to investigating approximate 
solutions of elliptic and parabolic second-order equations. Then, given 
an elliptic operator 

£ = dijP.Pj + hiDi 
with variable coefficients, a natural question arises as to whether it is 
possible to construct operators Lh such that they converge to £ and 
our assumptions are satisfied. 

This question has little to do with the dependence of a^- and bi on t 
and we assume that 

and a and b are bounded and continuous along with their first-order 
derivatives. 

It is not hard to see that under the symmetry assumption (S) the 
operators L h approximate £ with 

aij(x) = (1/2) £ qx(x)\i\j, (4.1) 
AeAi 

AGAi 

in the sense that L h ip — > £ip as h { for all smooth ip. 
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A few basic examples describing conditions on q\, which guarantee 
that our assumptions are satisfied, are given in Remarks 12.3} 16.4} and 
16.111 However, how this can be transformed into some conditions in 
terms of is not clear right away. 

By the way, it is shown in [1] that if C admits finite-difference ap- 
proximations constructed by contracting a fixed mesh and the approxi- 
mating operators satisfy the maximum principle, then they always have 
the form (12.21) with Ai = —Ai and q\ = q_\. This form is nonunique 
and the issue of choosing appropriate q\ and p\ arises. 

It is proved in [7J that, if the matrix a is uniformly nondegenerate, 
then there always exist Ai and q\ possessing property (S), such that 
(14.11) holds, q\ are as smooth as a is, and q\ > k > 0, where k is a 
constant. 

It is also proved in [7J that if all values of the matrix a lie in a closed 
convex polyhedron in the set of nonnegative matrices and a{x) has 
two bounded derivatives, then again there exist Aj and q\ possessing 
property (S), such that 04. ip holds, and ^/q~\ are Lipschitz continuous. 

In these two cases the issue of satisfying our assumptions reduces to 
representing b(x) appropriately. 

There is a way to do so, used quite often in probabilistic literature, 
by adding, if necessary, the set T = {ie 1 , ±e d } to Ai, where {e 1 } is 
the standard basis in IR d , defining 

(t± = (l/2)(|t| ± £)), and defining p\ = on the remaining part of 
r U Ai. There is a certain inconvenience in this approximation, which 
we discuss in the following example along with a way to avoid it by 
using different p\'s. 

Example 4.1. For d — 1 consider the operator 

Lu = bDu, 

where b = b(x) is a smooth function bounded along with its derivatives. 

If b changes sign, then, since in our setting xx — hp\ is required to 
be > 0, we have to take Ai consisting of at least two points {Ai, A3}. 
The most natural choice is Ai = {±1} and 

L h ip(x) = Ph<p(x) = b + (x)5 hjl <f(x) + b-(x)5 h - 1 ip(x). 

Observe that on smooth if we have Lhf — > bDip as h J. 0. 

Notice that if b changes sign, p±i = b± are Lipschitz continuous 
but need not be continuously differentiable unless we impose a severe 
restriction on the behavior of b near the points where it vanishes. Ac- 
tually, in this article the assumption that q\ and p\ are smooth can be 
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replaced with the assumption that they are Lipschitz continuous and 
then require (12. 3p be satisfied for each t almost everywhere with re- 
spect to x rather than for all x. However, in such case it is unrealistic 
to assume in Theorem 12.11 that u is continuously differentiable in x. 
Generally, u will be only Lipschitz continuous in x and estimate (12.41) 
will hold almost everywhere rather that everywhere in Qt- More seri- 
ous trouble occurs when we want to estimate higher order derivatives, 
which we will be concerned with in a subsequent article. Then we need 
Px to have higher order derivatives and this excludes many interesting 
cases. 

On the other hand, the reader may like to check that with the above 
p condition (I2.3P is satisfied (a.e.) with any 5 G (0, 1) as long as b is a 
decreasing function, which agrees well with the limit case of differential 
equations. 

One can construct a different approximation of bDip for which p\ are 
as smooth as b. Indeed, take a constant 9 such that |6| + 1 < 9 and set 

Pi = b + 9, p„i = 9. 

Then p± x > 1, again Pi^iV 9 + P-i^k-iP bDip on smooth (p, and 
p±i are as smooth as b. This method is somewhat close to adding 
an artificial diffusion. However, we add it only to the finite-difference 
approximation and not to the operator bD. 
Now consider the operator 

Lu = (l/2)aL> 2 + bDu, 

where we suppose that a(x) > and r := y/a and b are one time 
differentiable with derivatives uniformly continuous on IR. Again take 
Ai = {±1} an construct p\ as in Remark 16.61 and define q\ = a and 
r A = r = y/a. 

It is shown in Remark 16.71 that for h sufficiently small, Assumption 
12.31 holds (perhaps with different 5 and A'i), if 

14(r') 2 + 6'< (l-5)c + K ia . 

This condition describes what we need from £ in the one- dimensional 
case and it looks quite satisfactory. On the other hand, it is yet stronger 
than the common assumption 

|r'| 2 + b' < (1 - 5)c + K ia , 

which along with other standard assumptions guarantee that solutions 
of D t u = Lu — cu + f admit estimates of the first derivatives in x 
independent of the time interval. 
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5. Discussion of Assumption 12.31 

In a subsequent paper about higher order derivatives estimates we 
will impose the following assumption, which trivially implies Assump- 
tion 10 

Assumption 5.1. We have m > 1 and there exist a constant 5 G (0, 1] 
and an operator K, — JCh £ such that 

2m^(^)L° A T A ^ < (1 - S) Q(hv) 
AeA AeA 

+ K&ip) + 2(1 - S)cK( IM 2 ) (5.1) 

AeA 

on Ht for all smooth functions ip. 

In this section we are going to discuss Assumptions 12.31 and 15. II Here 
we suppose that only Assumptions 12.11 and 12.21 are satisfied. 

Remark 5.1. Condition (15 .ip involves a mixture of finite differences and 
derivatives. It is reasonable to split it into two parts, the combination 
of which turns out to imply (15. 1ft : For all smooth tp we have on Ht 
that 

2m (h<P)L°T x tp <(1-5)J2 Q(M 

AeAi AeAi 



+ KiQ((p) + 2(1 - 5)dC( l^l 2 ), (5.2) 

AeAi 

that is (15.11) holds with tq = and 

d 

2mJ2(D i <p)L° i cp < K i Yl Q ^¥) 

i=l AeAi 

+ (l/2)5c\D ¥ >\ 2 + KM I^H, (5-3) 

AeAi 

where and below by /C we denote generic operators (perhaps, depending 
on h) of class & and 

L° l¥ > = Y,{h- l D iqx + D iVx )8 x y. 

AeAi 

To show that (15.21) combined with (15. 3p imply (15.11) if we choose small 
r > appropriately, observe that in terms of Si. = TqDi equation (15.31) 
means that 

2m (M^aV < K ir % G(M 
AeA 2 AeAi 
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+ (l/2)*c l^l 2 + Kir 2 M \*M 2 )- (5-4) 
agA 2 AeAi 

By choosing r so that K\Tq < 5/2 and K\Tq < (1/2)Scq, slightly re- 
defining K, to absorb the second term on the right in (15.41) , and summing 
up (15. 4p and (15. 2p we come to an inequality which is even somewhat 
stronger than (15.11) if 5 there is replaced with 5/2, which is irrelevant. 

Remark 5.2. Assume that the symmetry condition (S) holds, q\ > 0, 
Px > 0, and r x := ^/q\ are Lipschitz continuous in x with a constant 
independent of t. Then it turns out that condition (15.31) is satisfied 
with any 5 G (0, 1), T\ = 1, and appropriate K\ and unit K. 

To show this observe that, for any unit £, \q\{p)\ < Nr x with N 
being the doubled Lipschitz constant of r\. In particular, by Holder's 
inequality 

( £ q m A^) 2 <N{J2 V~%\^M) 2 < N E (5.5) 

/xeAi ^eAi /jgAi 

which allows us to make obvious changes in the estimates of 1% in 
the end of Section El one of the changes being that now we can allow 
Q(5\ip) to enter the estimates with as large constant as we wish. 

In the following remark we discuss an estimate which was crucial 
in the nonlinear setting for establishing a rate of convergence of dif- 
ference approximations to the true solutions of Bellman's equations in 
cylindrical domains (see [3]). We will see how using different T\ can 
help. 

Remark 5.3. Consider the situation when the coefficients q\, p\ and c, 
the free term /, and the terminal data g also depend on a parameter 
y e R: 

q x = qx(t,z), p x = p\(t,z), c = c(t,z), f = f(t,z), g = g(z), 

where z = (x,y) G M. d+l . Assume that these functions and their first 
derivatives in z are bounded on H' T = [0, T] x M. d+1 and continuous 
in z. Assume that c > cq for all values of the arguments. Suppose 
that (15.21) holds on for m — 1, T\ = 1, any smooth <p(x), and any 
value of the parameter y with K, perhaps depending on y (as well as h 
and t). Assume also that the symmetry condition (S) holds, q\ > 0, 
P\ > 0, and r\ = ^/q~\ are Lipschitz continuous in z with a constant 
independent of t. 

Finally, suppose that in H' T we are given a bounded function u(t, x) = 
u(t, x, y) which satisfies (12 .ip in Q T for each value of y. Of course, now 
in (12.11) we write z in place of x. 
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Take an e > and set 

T*u(t,x,y) = u(t,x,y + s), 5* = e' 1 ^ - 1). 
We claim that in [0, T] x Q x R it holds that 

\S»u\<N( sup + sup + V |£ A u|)), (5.6) 

H T xR (S'Q T )xR ^ i 

if £ G (0, r/i], where TV and r G (0, 1] depend only on 5, c , i^i, the 
number of elements in Ai, |Ai|, and the Lipschitz constants of c, r A , 
and p\ with respect to z — (x, y) . 

To prove the claim observe that, although equation ( 12. ip can be 
considered in as an equation with parameter y, we will treat it as 
an equation in Q' T = Qt x M. Then we denote by Ao the positive vector 
on the y-axis having the length e/h, take a r > to be specified later 
and introduce 

A' x = Ai U {A }, q Xo = p\ = 0, r Ao = rh/e, r x = 1, A G A x . 

We now check that Assumption 12.31 is satisfied for the new objects 
with m — 1, To = and A[ in place of A. Owing to the assumption 
that (I5.2p holds, we immediately see that the left-hand side of (12.31) for 
the new objects is less than 

(1 - 6) J2 2(M + K x Q(ip) + 2(1 - 8)cK,( ^ \5 xV \ 2 ) + 7, 
AeAi AeAi 
where I = I\ + I2, 

Since for smooth ijj, 

5 Xo ip = T5 y £ ip, \5 Xo ip\ < rsup \dip/dy\, 

we have that 

\I 2 \<NtJC(J2\5^\ 2 ). 

(ueA'j 

Upon observing the following general properties of finite-differences: 
h5 x 5, = (Tp - 1)5 X , A,T X = -5,5., + 5 X 5, + 5 X 5^ (5.7) 
and combining them with the estimate 

r- l \Sx Q,\ = = |2r M 5 £ % + e{5 y £ r,) 2 \ < + h), 

and ( 15.51) . one easily shows that 

|/i|<iVr(^Q(^) + /C(^|^| 2 )). 
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Since e < rh, we have that \5\ip\ < \5xv\ on A[ and the above estimates 
show how to choose r G (0, 1] in order for Assumption 12. 31 to be satisfied 
indeed. 

By applying Theorem 12.11 (and Remark 12.41) we finish proving our 
claim. The only point which is perhaps worth noting is that by The- 
orem [2J] the constant N from (15.61) also depends on \A[\. However, 

|a;|2 = | Al |2 + r 2. 

6. Discussion of Assumptions 12.31 and 15.11 in case that 

t x = 1,A g Ai 

Here we suppose that only Assumptions 12.11 and 12.21 are satisfied. 
Everywhere below we set T\ = 1 for all A G Ai. 

Remark 6.1. Suppose that the symmetry assumption (S) is satisfied. 
Then the operators can be regarded as finite-difference approxi- 
mations of £ (see (12.81) ) in the sense that for any smooth ip we have 
Lhf — > dp as h I 0. If we are only interested in this property, then we 
can always assume that p\ > K\. 

Indeed, if we do not have this inequality, then we take a sufficiently 
large constant K 2 (independent of h), redefine p\ as p\ + K 2 . This will 
not violate the convergence Lh<p — > Cp since 

AeAi AeAi 

if tp is smooth. 

The following lemma is often used below and in the continuation of 
the present paper. 

Lemma 6.1. Let be a nonnegative function on ~R d for each /i G Ai 
and A from a finite set of indices A'. Assume that 

sup a v < C 

for some function C onM. d . Then there is a JC G A such that 

MeAi,AGA' AeA' 

on ~R d for any bounded Borel function f = fx given on M. d . 
Proof. Using 

h\8,f x f = ((T„ - 1) fx) 2 < 2{T,fxf + 2fl = 2(T, + 1)/ A 2 , 
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we can estimate from above the left-hand side of ( 16.11) by 

2 E C ^ + !) E fl 

/xeAi AeA' 
where := sup AgA / a^. Hence we get (16.1 p with 1C e & defined by 

w) = h E w + !)/■ 

/aeAi 

The lemma is proved. 

Remark 6.2. One can give sufficient conditions for (15.21) to hold without 
involving test functions if, which makes them more "explicit" and in 
combination with Remark 15.21 covers many situations when Assumption 
15.11 is relatively easy to check. One set of these "explicit" conditions is 
given in this remark. Here we also show why the operators of class ^ 
are useful and how the presence of hp\ in xx entering the operator Q 
on the right of (15. 2p may help. 

Suppose that Ai = — Ai and q x = > and set r A = Take 
a 5 G (0, 1/4) and assume that on H T there are functions r AA1 , p\^ > 0, 
A, \i G Ai such that 

h 2 {5 x r,) 2 < 5( X , + Xx) + h 2 r 2 Xfl , E SU P r l ^ 25c ' ( 6 - 2 ) 

^ AeAl 

h 2 \5 x p^\ < 5 2 (xn + Xx) + Sh 2 p x ^, E SU P P>* - ( 6 - 3 ) 

T7 AeAi 
/ueAi 

By virtue of Remark 16.11 if Lh are used for approximating C, we can 
change these operators and have (16.21) and (16.31) satisfied with r Xfi = 
PX[i = for sufficiently small h, provided that the Lipschitz constants 
in x of r x and p\ are bounded with respect to t. 
For a function £\ given on Ai let us write 



l£l 2 =E le 



x\ 

xeAi 



Then it turns out that condition (15.21) is satisfied if on for all func- 
tions £a w e have 

10m 2 (l - A5)~ 1 Ji + 2m 2 (l - 4<TT 1 J 2 
+25m 2 E ex\5xP,\+2m E ZxUhp» + (6xr») 2 ) 

< (2 - 85)c|£| 2 + Kj E ^aXa + Sh~ 2 E Xa|6 + C-a| 2 , (6.4) 

AeAi AgAi 
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where 

To prove this, use formulas (15. 7\i . Also drop the summation sign 
having repeated indices in Ai to see that 

2(5x(p)Q\T x p = (S x ^)(Sxq^)(Sx^ + <5 A <5_^ + A^Jy? 

= 2£ A (<5 A <^)<5 A <5 M y9 - Cxihq^^-^ = hi + I12 + hi + I22, 
where £ A = 8\<£ and 

Ai = ^x(hr^)r^5 x S^ip, I 12 = 2h£ x (5 x r^) 2 5 x 5 l ,ip, 

hi = -Ztixihr^r^-^nV, I22 = h^ x (hr^) 2 ^^ = 2£ A (<5 A r M ) 2 £ A1 . 
Before starting to estimate 1^, we note that as h J, the terms In, I 22 , 
and /12 disappear if if is twice continuously differentiable due to the 
symmetry of Ai. In that case there is no need to estimate them. For 
fixed h they are present and estimating them is only possible under 
stronger assumptions than in the case of partial differential equations. 

Now notice that by Young's inequality 

ml n < (1/2)(1 - 45)/ + 8m 2 (l - 45)~ l J u 

where 

/ = £ Q(s xV ) = £ x.lWM 2 . 

AeAi A,^t6Ai 

Similarly, 

mhi < (1/2)(1 - 45)/ + 2m 2 (l - A5)~ l J 2 . 
Next, owing to (16.21) 

m/12 = 2m(£ A <5 A r M ) (h^S^^S^^) 

< 2m 2 J x + (l/2)/i 2 £ (5 A r M ) 2 (5 A 5^) 2 
A^eAi 

< 2m 2 Ji + (5/ + (l/2)/z 2 ]T r 2 M ((5 A <5^) 2 , 

A,^eAi 

where the last term by virtue of Lemma 16.11 is estimated by 

4<5c/C(]T((5 A ^ 2 ). 

AeAi 

By collecting the above estimates we obtain 

2m(5 x <p)Q x T x ip < (1 - 35)/ + 10m 2 (l - 4<5)- 1 J 1 

+ 2m 2 (l - AS)' 1 J 2 + 2m{5 x r ll fU„ + 4c5c/C( ]T ((5 A <^) 2 ). (6.5) 

AeAi 
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Next, 

2(6 x <p)P x T x <p = 2Uhp^ + 2ht x (6 xPli )6 x 6 li <p, 
where the last term is maj orated by 

25m £aI<W + (lmm^S-W £ \5 xP ,\{5 x 5^f . 

A,/iGAi A,/iGAi 

We use assumption (16. 3p and proceed as while estimating 1\2- Then we 
see that 

2m(S X ip)P x T X ip < 2m£x(8xPr)£p 
+ 25m 2 £\ S *P»\ + 51 + 25c/C(^(5 A ^) 2 ). (6.6) 

A,/igAi AgAi 

Finally, upon combining (16.61) with (16. 5p we obtain 
2m(8 X (p)L° x T x ip < (1 - 25)/ + 10m 2 (l - Ad)' 1 J x + 2m 2 (l - 45)" 1 J 2 

+25m 2 e x \hp»\+6ScJC(J2( 6 ^ 

A,/iGAi AgAi 

+ 2me A ((5 A r M ) 2 + (5 A ^))e M . (6.7) 

We use the fact that h~ l (5 x + 5- X ) = — 5 X 5~ X , use assumption (16. 4(1 . 
and take into account that for any /C' G .ft 

(2 - 85)c IM* + 65c/C'( ]T |M 2 ) = (2 - 25)c/C( £ |5 A ^| 2 ) 

AgAi AgAi AgAi 

with an appropriate K G A. Then we estimate the right-hand side of 
(1-5)1 + KiQ{ip) + (2 - 25)cK{ £ (c^) 2 ), 

AGAx 

and we see that (I5.2p is satisfied indeed. 

Remark 6.3. It is easy to see that if (16.21) and (16. 3p hold and we assume 
that inequality (16.41) is satisfied with an additional term ty(t,x,£) on 
its left-hand side for some function \P of t, x and £ = (£a)a6Ad then 
inequality (15. 2p holds with the additional term ^f(t, x, (5a<£>)agAi) on its 
left-hand side. 

Remark 6.4. Assume (S), assume that q x > 0, p x > and let m > 2. 
Then it turns out that Assumption 12.31 is satisfied for 5 = 1/10 and 
appropriate ro > or r = if Co is sufficiently large (independently of 
h). 

Indeed, it is well known that the Lipschitz constant of the square 
root of a nonnegative twice continuously differentiable function w(x) is 
controlled by the supremums of its second order derivatives. Therefore, 
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Remark I6T21 (where we take r AA1 = <5 A r M andpA^t = 10|<5aP/J) immediately 
implies that condition (I5.2p is satisfied if c is large enough. 

That condition (15.31) is satisfied follows from Remark 15.21 and again 
from the fact that for twice continuously differentiable w on M. d we have 

l-Du'l 2 < 2w sup 

xGM d , |C|=1 

At this point we do not even need large cq. Referring to Remark 15. II 
we obtain what we have claimed. 

Actually, above in this remark we used that m > 2 only to guarantee 
that r\ are Lipschitz continuous in x with a constant N' independent 
of t. If we just assumed this last property, then our argument about 
(15.21) would become even shorter. In addition, what was said about 
(1531) is still valid. 

Remark 6.5. Under the symmetry assumption (S) and the assumption 
that q\ > and p\ > one can give a rougher condition without using 
£a and implying (16.41) with m = 1 and sufficiently small 5. Then ( 15.21) 
will be satisfied as long as conditions (16.21) and (16.31) are. 

By the way, recall that, for all small h, one can always satisfy condi- 
tions (16.21) and (16. 3p on the account of modifying if necessary p\ if the 
Lipschitz constants of r\ in x are bounded in t (see Remark 16. 2p . 

By the inequality 

Y V\V»a\» < Y Y \ ax » + a ^ 

A.^tgAi AeAi M^Ai 

we have 

E ^ E ^ E + ^ A '' 

A,^eAi AeAi /iGAi 

J2= y E owow ^ 2 Y & E i E oww- 

\,u,eAi AeAi i^eAi /xeAi 

A simple argument based on the above estimates and continuity 
shows that (16. 4p holds with m = 1 and a small 5 > if for any A G Ai 

io Y(^,) 2 +^Y\J2^ r ^\ 

AteAi i/eAi ^eAi 

+ 2 Y l<ton + ^ + (5a^) 2 + (<Va) 2 | < c + K\q\. (6.8) 

MeAi 

Condition (16. 8p basically means that if for a A G A x at some point in 
H T the value q x is small, then either <5a?V, 5^t a, 5 m £>a, an d 5\p^ should 
be small or c be large at this point. As the point varies, the dominating 
terms may change roles. 
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Remark 6.6. There are cases when it is preferable to keep the last term 
on the left in (16.41) as is. 

To see a reason for that, let m = d = 1, take a constant 6 > 1 and 
define Ai = {±1}, 

g A = 0, p 1 {t,x) = {l/2)b{x) + 9, p_i = -{l/2)b{x) + 9, 

where b(x) is a decreasing function with bounded derivative such that 
|6| < 1. Observe that for any 9 we have L° h (p — > bip' as h J. if ip is 
smooth. 

Now notice that condition (16.21) is trivially satisfied since r\ = 0. 
Condition (16.31) is satisfied with any 5 > and p\^ — if h < 1 and 
> 100 sup |6'| + 1 since 

Xx = hp x >h(9-\b\)>h(9-l) 

The left-hand side of (16. 4p is one half of 

AS £ &\5M + (6 - ^i) 2 [^i& - + (£ 2 - + 5-ib). 

i=±i 

Here the middle term is nonpositive since b is decreasing. Also the last 
term is maj orated by 

Furthermore, concerning the right-hand side of H6.4I) observe that 

A A A 

It follows easily that, no matter how small Co is, for sufficiently small 
5 and large 9 condition ( 16.41) and, by Remark 16.21 condition ( 15.2ft are 
satisfied. This along with the almost obvious fact that (15. 3p holds 
shows that Assumption 12.31 is satisfied as well. 

To finish the remark notice that if we tried to check condition (16.81) 
we would fail to do that for small h no matter how large 9 is unless c 
is large enough. 

Remark 6.7. We continue the analysis of the one- dimensional situation 
started in Remark 16.61 So, we assume that d — m — 1 and we have in 
mind approximating an operator C<p(x) = (l/2)a(x)<p"(x) + b(x)<p'(x). 
As in Remark l6.6l we assume that |6| < 1 and, in addition, assume that 
a > and r := yfa and b are one time different iable with derivatives 
uniformly continuous on E. 

As in Remark 16.61 take Ai = {±1} and define p±i for an appropriate 
9. Then also set r M = r. Now observe that both parts of (16.41) are 



28 I. GYONGY AND N. KRYLOV 

order-two homogeneous functions of £\. Therefore, it suffices to check 
(16.41) assuming that 

E £ = I- (6-9) 

AeAi 

Then, due to the assumption that r' and b' are uniformly continuous 
on M, it is not hard to see that 

J 1 ~2(r') 2 , J 2 ~2(r') 2 (£i-£-i) 2 , 

where by a ~ (3 we mean that for small ft. the difference |a — /3\ can 
be absorbed into c > Co > with as small coefficient as we wish. 
Therefore, upon recalling the estimates from Remark 16. 6} we see that 
condition (16.41) is satisfied for all small h if 

20(1 - 45)~ 1 (r') 2 + 4(1 - 45)- 1 2(r') 2 (£i - £-i) 2 + 2(r') 2 (£i + £-i) 2 

+ 25|6 / | + (6-e-i) 2 & / < (2-95)c+K 1 a + 5h'\9-l)^ 1 + ^ 1 ) 2 . (6.10) 
On the account of (assumption (16.91) and) the presence of h~ l (6 — 
l)(£i + £-i) 2 on the right of (I6.10p . it suffices to check (16.101) for small 
h assuming that the inequality |£i + < h 1 ^ 2 holds. It follows that 
(IBTTOj) holds for small h if it holds with (2 - 105)c in place of (2 - 95)c 
but only for £\ satisfying £i = — In that case (£i — £-i) 2 = 4^ = 2 
and (I6TT0|) holds if 

28(1 - 4<5)" 1 (r') 2 + 26\b'\ + 2b' < (2 - 10<J)c + K x a. 

Since we would be satisfied if ( 16. 4p held with at least one 5 > 0, we see 
that, under the assumption of the present remark, (16. 4p (perhaps with 
different 5 and Ki) is indeed satisfied for small h if 

14(r') 2 + 6' < (l-5)c + K ia . 

Remark 6.8. There are multi-dimensional analogs of the situation in 
Remark [6761 For instance, let U(x) be a concave function with bounded 
derivatives and assume that 

Pfl = (DU,n) + 6, 

where 6 > max^gAj \fi\ sup \DU\, so that > 0. Then for small h we 
have 

E &£Apm~ E U^ 2 ^,A) = (d 2 ^^)<o, 

A,^eAi A,/^eAi 
where 77 = ^2 
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Remark 6.9. Recall that 5\ = 5h,x and assume that ( 15. 3p holds for all 
small h > with /C perhaps depending on h. Then it turns out that 
for all t £ [0, oo] 

^^A<7a(£, a;) is independent of x. (6-11) 

AeAi 

To show this observe that, since the values of the first derivatives of 
(f at a fixed point have nothing to do with the increments of (f, ( 15.31) 
is equivalent to saying that 
d 

j2(l^) 2 <n( xa|^i 2 +/c(^i^i 2 )), 

i=l X,ueAi AgAi 

where the constant N can be easily computed given supc and K\. It 
follows that 

d 

J>W < N( £ Xx\Sx6M 2 + sup l^l 2 )- (6.12) 

i=l A,i/GAi Rd AeAi 

Now multiply (16.121) by h 2 and let h [ 0. Then we obtain 

d 

V [(D^fijjVA^tt,^] 2 < 7Vsup|^|. 
i=i /igA Kd 

which leads to the conclusion that for smooth ip, % = 1, ...,<i, and all 
A £ A 

mga 

This is equivalent to saying that (16.111) holds. 

Remark 6.10. Additionally to Assumptions 12. l l and 12.21 suppose that 

A + z/£"Ai VA,i/eAi. (6.13) 

It turns out that in this case (I5.3P is satisfied for all small h only if for 
any A G Ai 

(i) either —A £ Ai and 5a (t, = q~\(t, x) +r\(t) for a function r^(t) 
independent of x, 

(ii) or — A £" Ai and q\ is independent of x. 
In particular, 

g^ = (l/2) Yl M^xf E- :=0). (6.14) 

AGAin(-Ai) 

We may concentrate on proving our claim assuming that q\ is inde- 
pendent of t. As we have pointed out in Remark 16.91 condition (15.31) 
implies ( I6.12p . We write the latter at x — 0, substitute ip(x/h) in place 
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of ip and let h I 0. Then by just comparing the powers of h in different 
terms we obtain that for a constant N and any i = 1, d and tp we 
have 

E (A?a(0)MA) < N( E |5 liA ^(0)| 2 ) 1/2 . 

AeAi XmgAj 

We see that the linear function (of ip) on the left provides a supporting 
plane at the origin for the convex function on the right. Consequently, 
there are some constants q\^ such that for all ip 

E(A<?a(o)MA) = E 9a,AaVv(o) 

AeAi A,i/eAi 

= e <^x A + v ) + ?(°) - ^( A ) - ( 6 - 15 ) 

A,^eAi 

Without losing generality we may assume that 

Qx,u = Qv,x 

and split the sum on the right in (16 . 1 5[) into two parts: the first part 
with the summation over A, v such that A + v = and the second part 
for X + u 0. According to assumption (I6.13P the terms p{\ + v) in the 
second part do not appear elsewhere in (16.151) . It follows that q\ M = 
if A + v ^ 0, so that 

E (Afc(O)MA) = E mW) - v(A) - v^(-a)]. 

AeAi AeAi 

Here the expression on the right is symmetric with respect to the trans- 
formation <p(x) — > <p(—x). Thus, 

E(A<?a(0)V(A) = £(A?a(0)M-A). 

AeAi AeAi 
We obtained this relation at the origin. Similarly, for any x 

E (A?a(x)MA) = E (A<?a(x)M-A). (6.16) 

AeAi AeAi 

Fix a Ao G Ai an take a <p which is 1 as A = Ao and zero otherwise. 
Then (1636) shows that 

(i) either — A £ Aj and then Dq\ (x) = Dq^\ (x), 

(ii) or — A £" Ai and then Dq Xa (x) = 0. 
This proves our claim. 

In Remark 16.101 we saw that (I6.13P along with (15. 3p lead to the 
symmetry of the operator Q u expressed by (I6.14p . However, if // + A £ 
A x for some /i, A £ A l5 the symmetry of Q u may not occur. 
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Example 6.1. Let d = 1 and Ai = {—3,-1,1,2}. Take a smooth 
function f(x) such that 1 < / < 2 and set 

g_ 3 = 1 q-i = q2 = 3 - /, gi = /. 

Then 

]T q x X = -3 - (3 - /) + / + 2(3 - /) = 0. 

AeAi 

Furthermore, 

h- 1 Yl = h ' 2 f'(- T -i + 1 + T i ~ T 2)V = f'Bp, 

AeAi 

where 

R:=h~ 2 (T 2 -l)(T_ l -l) = 5 2 5^ l . 
It follows that (15. 3p is satisfied. Also observe that for \i G Ai we have 

AeAi 

and 

AeAx 

This and the fact that q\ > 1 easily imply that condition ( 15. 2ft is also 
satisfied with appropriate constants and operator /C in case there are 
terms also with p\ in L and either px > or h is sufficiently small so 
that xa > 1/2. 

Remark 6.11. Note that the argument in the above example shows that 
(always under Assumptions 12.11 and 12.21) an operator 

L h = hT 1 Y <?a5a - c 
AeAi 

satisfies Assumption 12.31 if q\ > n > for a constant k > and the 
equality 

h- 1 J2(D l q x )5 x = <fcvA^ i = l,2,...,d (6.17) 

AeAi A,/ieAi 

holds with some bounded coefficients q^. Therefore it would be useful 
to find simple conditions, i.e., which can be easily verified, for the 
characterization of Ai and q\ satisfying (16.171) . In this direction we 
have the following condition and conjecture about a criterion for (16.171) 
to hold. 

We call a function ip on Ai U {0} linear if ip(v) = ip(X) + 
whenever v, A, fi G Ai U {0} and v = n + A. 
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Conjecture. Equation (I6.17P holds with some if and only if 

AgAi 

holds for any if that is linear on Ai. 

Notice that (I6.17P and the condition of the conjecture are satisfied, 
for example, when property (S) holds, or if Ai is the union of disjoint 
triplets {A, /i, A + such that —Dq\ + ^ = Dq x = Dq^ for each of them. 

Remark 6.12. Condition (15. ip is "fool proof in two ways related to 
changes of variables. For simplicity we only concentrate on the case 
that ro = and T\ — 1 for A G Ai. First, one can try changing 
the time variable by introducing the new function v(t, x) = u(k~H, x), 
where k > is a constant. This amounts to dividing the coefficients of 
(12.11) and / by K and accordingly changing time t — > n~ 1 t. However, as 
is easy to see this will not affect condition (15. ip and, for that matter, 
the value of F\ in Theorem 12.11 either. 

The second way is to try to relax condition (15.11) by changing the 
space variable. Introduce 

h = Kh, S : ip — > Si/}(t, x) = i/)(t, k^x), 
qx = n 2 Sq x , px = KSpx, 

= Y v x6h > x + Y p^h,\, 

AgAi AgAi 

£}U = Y ( 5 h,^) 5 h,X + Y ( S hvP><) 5 h,\- 
AgAi AgAi 

One easily checks that 

8- hjX = K~ 1 S5 h) xS~\ AeA, L\ = SL h S-\ (6.18) 
Owing to (I6.18p . if u satisfies 00 (12.11) . then u = Su satisfies 

— = L- h u-cu + f, 

where c = Sc, f = Sf. Furthermore, = k _1 1 SX/ 1>/j S'~ 1 and if (15.11) 
is satisfied, then for ip = S^ip 

AgAi AgAi 

<(1-5)k~ 2 S (qx + h Px )\5 htX 5 ht ^\ 2 

i/GA.AeAi 

+K 1 k~ 2 S (<?a + hpx)\6 h ,xi>\ 2 + 2(1 - 5)K~ 2 S[c1C h { ]T \5 h ^\ 2 )] 

AgAi AgA 
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f£A,A(=Ai 

+K.K- 2 J2(q x + ^ A )|<Mf + 2(1 - VcSKhS- 1 ^ \k,M 2 ), 
AeAi AeA 
where SJChS~ l G .ft. We see that this change of coordinates did not 
produce any effect on (15. ip apart from changing K\ and JCh, which is 
irrelevant. 
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